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4 Abstract. Neutrinoless double beta decay is a lepton number violating process whose 

5 observation would also establish that neutrinos are their own anti-particles. There are many 

6 experimental efforts with a variety of techniques. Some (EXO, Kamland-Zen, GERDA phase I 

7 and CANDLES) started take data in 2011 and EXO has reported the first measurement of the 

8 half life for the double beta decay with two neutrinos of 136 Xe. The sensitivities of the different 

9 proposals are reviewed. 
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O , io 1. Introduction 

For many isotopes like 76 Ge j3 decay is energetically forbidden, but double beta decay {2vf50) 
is allowed 



{A, Z) -> (A, Z + 2) + 2e~ + 2z7 e (1) 

11 This was suggested very early pQ and - following the idea of Major ana that neutrinos could 

12 be their own anti-particle [2j - also the possibility of neutrinoless double beta decay was 

13 anticipated shortly afterwards [5] (for a review see [HE]). The latter case is very interesting 

14 since lepton number is violated and it would establish that the neutrino is its own anti-particle. 

15 The experimental signature in this case is a line at the Qpp value of the decay if the sum of the 
O i6 electron energies is plotted. 

17 Searches for double beta decay date back to the beginning of nuclear physics and nowadays 

is more than a dozen large scale experimental programs are suggested. These programs are 

19 compared in this article and also the status of theoretical matrix element calculations is 

^ 20 discussed. For general reviews the reader is refered to the literature, e.g. [HI El [8]. 

21 There are also other related processes like double positron decay or double electron capture 

22 processes. While 0i//3/3 is already a suppressed process, the other decays are expected to be even 

23 rarer unless there is some resonance enhancement [9t ll0l[TT| ll2j. In this article only decay 
searches are discussed. 
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25 2. Motivation 

26 The observation of neutrino oscillation establishes that these particles have mass [13]. Since 

27 neutrinos have no electric charge, there is no known symmetry which forbids additional terms 
in the effective Lagrangian beside the standard Dirac mass term mo [Zj: 
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L Yu k = m D v L v R + ^m L v L (v L ) c + -m R {y R ) c v R + h.c. (2) 
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{pZi {vr) c ) " u " nJ )+h.c. (3) 



29 The subscript L stands for the left-handed chiral field vl = |(1 — 75) v and R for the right- 

30 handed projection |(1 + 75) v. The superscript C denotes charge conjugation. The jur term 

31 describes therefore an incoming neutrino and an outgoing anti-neutrino, i.e. lepton number is 

32 violated by 2 units. The eigen states of the mass matrix are of the form (y + v c ). Consequently 

33 neutrinos are expected to be - in general - their own anti-particles, i.e. Majorana particles. 

What is the best experimental approach to establish that our known neutrinos are Majorana 
particles? Neutrinos (or anti-neutrinos) are produced in charged weak current reactions and 
- depending on the charge of the associated lepton - only one chiral projection couples. For 
example in (3 decay n — > pe~ v e ,Rt a right-handed anti-neutrino couples: 

1 3 

V e ,R = "e^i 1 + 75) = 2j U ei(^i,h=+l + -^^i,h=-l) ( 4 ) 
~ i=l 

34 Here, U is the PMNS mixing matrix, vi are the mass eigen states with mass eigen values m,;, E 

35 is the neutrino energy and h stands for the helicity of the anti-neutrino. 

For a Dirac particle these anti- neutrinos can only undergo detection reactions like p v e ,R 
ne + . If, on the other hand, neutrinos are Majorana particles, then the v^h = —\ component can 
undergo the reaction v e l n — >■ pe~ with 

1 3 m 

Ve,L = - (1 - lb)v e = U ei{^i,h=-1 + ^H,h=+l) (5) 
i=l 

36 The rate of this reaction^] is however suppressed by the factor (rrii/E) 2 which is e.g. 10~ 14 for a 

37 neutrino mass of 0.1 eV and a neutrino energy of 1 MeV. Thus solar neutrino experiments for 

38 example will not be able to establish the nature of neutrinos. 

39 The alternative is the search for 0v(3/3 where the neutrino only enters as a propagator 

40 ~ mpfi/q 2 = Y^i C^ei ' m i/Q 2 - The coupling strength is called the effective Majorana mass. 

41 Since one mole contains a large number of nuclei, the factor (rrii/E) 2 is compensated. For 35 

42 isotopes double beta decay is the only possible decay mode. The Standard Model allowed decay 

43 with two emitted neutrinos (2f/3/3) has been observed for 11 isotopes with half lives between 

44 7 • 10 18 y and 2 • 10 21 y O [TS]. 

45 Part of the Heidelberg-Moscow experiment claims to have observed 0v(3(3 of Ge with 

46 mpp ~ 0.2 — 0.6 eV [16J. Clearly this needs independent confirmation which poses another 

47 motivation for the experimental efforts. 

48 3. Experimental sensitivity 

An experiment will observe some background events Abkg which 
detector mass M - is given by 

Abkg = Af ■ t ■ B ■ AE 

and possibly signal events 

A sig = In 2 • N A ■ e ■ r, ■ M ■ t/(A • T^ 2 ). (7) 

49 Here t is the measurement time, B the so called background index given typical in 

50 cnts/(keV-kg-y), AE is the width of the search window which depends on the experimental 

51 energy resolution, Na is the Avogadro constant, e the signal detection efficiency, r\ the mass 

52 fraction of the 0v (3/3 isotope, A the molar mass of this isotope and T®Y 2 its half life. 



if this number scales by the 

(6) 



1 Since the charge of the outgoing lepton is the same as in the production process, U and not U* enters here. 



Table 1. List of most interesting Ov(3(3 isotopes. Half lives are taken from |14| [15] while all 
other numbers are from [7J. 



isotope 


qOu 

[^] 


[keV] 


nat. abund. 

[%] 


rn2l/ 

1 l/2 

[10 20 y] 


experiments 


48 Ca 


6.3 


4273.7 


0.187 


0.44 


CANDLES 


76 Ge 


0.63 


2039.1 


7.8 


15 


GERDA, Majorana Demonstrator 


82g e 


2.7 


2995.5 


9.2 


0.92 


SuperNEMO, Lucifer 


100 Mo 


4.4 


3035.0 


9.6 


0.07 


MOON, AMoRe 


116 Cd 


4.6 


2809 


7.6 


0.29 


Cobra 


130 Te 


4.1 


2530.3 


34.5 


9.1 


CUORE 


136 Xe 


4.3 


2461.9 


8.9 


21 


EXO, Kamland-Zen, NEXT, XMASS 


150 Nd 


19.2 


3367.3 


5.6 


0.08 


SNO+, DCBA/MTD 



If Abkg < 1 the experimental sensitivity scales with M • t while for Abkg >> 1 the e.g. 90% 
C.L. limit on the half life (assuming there is no signal) is given by 

If systematic errors become important e.g. if the energy resolution is not well known or the 
assumption of the background shape is not correct, then the sensitivity is reduced. 

4. Theoretical considerations 

The half life for 0uf3f3 is given by [7] 

2 

[T^]- 1 = GouiQpp, Z)\M 0u \ 2 ^f (9) 

Here Gq v is the calculable phase space factor (Tab. [I]), m e is the electron mass and Mq v is the 
nuclear matrix element whose calculation is difficult and can only be done using approximations. 
For a review see for example [U [8] . 

While the observation of Oz^/3/3 would manifest lepton number violation and the neutrino's 
Majorana nature, the underlying physics can only be disclosed if the observed T®J 2 for different 
isotopes and possibly other variables like the angle between the emitted electrons is compared 
to theory. Consequently, there is a large interest in nuclear matrix element calculations 
and substantial progress has been made during the last years. Traditionally, nuclear shell 
model (NSM) and quasi particle random phase approximation (QRPA) calculations have been 
performed. Recently new approaches like the interacting boson model (IBM), the generating 
coordinate model (GCM) and the projected Hartree-Fock-Bogoliubov (pHFB) method have been 
applied. A discussion of these calculations is given in |27j . 

The results of these calculations are shown in Fig. [TJ The following statements can be made 
concerning the status: 

• There is no large variation for the NME between the different isotopes. This might be due 
to the fact that only neighboring neutrons in a nucleus contribute to the decay [17\ [20] . 

• For the NSM, all values are systematically lower than for other methods. Possible reasons 
for this effect are discussed in the literature |1 7|, [2"8~] . 

• The differences between the QRPA calculations of different groups are now quite small. 
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Figure 1. Left: Different calculations for nuclear matrix elements for decay for light 

neutrino exchange. NSM = Nuclear shell model [HI [18], SRQRPA = self-consistent renormalized 
quasi-particle random phase approximation |19 [ 120 1 f2Tj (matrix elements are scaled by 1.14 to 
compensate for different phase space factors), pnQRPA = proton- neutron quasi particle random 
phase approximation [22J, GCM = generating coordinate method [23J, IBM = interacting boson 
model |24[ [25] (matrix elements are scaled by 1.18 to estimate the effect if the UCOM short 
range correlation instead of the Jastrow type would have been used |8j), pHBF= projected 
Hartree-Fock-Bogoliubov model |26| . Right: ratio of expected events per kg target mass 

for the different models normalized to 130 Te. 



• For a given isotope the calculations spread by typically a factor of 2, i.e. a factor of 4 for 

1 l/2- 

• The role of short range correlations has been studied and the UCOM correction has 
emerged as standard |29|. Alternatively, a self consistent implementation was first applied 
to SRQRPA [20 1 and later to other methods |18|. [26] and resulted in small changes. 

• Experimental input can have a large shift of the result. For example charge exchange 
reaction measurements of 150 Nd( 3 He,t) and 150 Sm(t, 3 He) |30] result in a quenching factor 
of 0.75 for the qa coupling and hence a reduction of the matrix element by 25% for 150 Nd 
|21j . In this calculation, deformation was treated for the first time in a QRPA calculation. 
For 76 Ge and 76 Se, the proton and neutron valence orbital occupancies have been measured 
[31] [32] . If the models are adjusted to reproduce these values, the NSM result increases by 
15% p.8J while the QRPA results are reduced by about 20% [33l|3l]. Hence the difference 
between NSM and QRPA becomes half as large. 

The calculations are performed for the standard light neutrino exchange but results for other 
mechanisms like SUSY particle exchange are also available |35[ [19] . 

In order to see whether some isotopes are better suited for 0i/(3(3 decay searches from a 
theoretical point of view, the number of expected decays per isotope mass can be compared. 
This value includes the phase space factor, the matrix element and the mass number A. For 
the comparison it is sufficient to look at the ratio of decay rates and in this case, some of the 
systematic effects of the matrix element calculations cancel since there are typically correlations 



Table 2. Selection of Ois/3/3 experiments. 



experiment 


isotope 


mass [kg] 


method 


location 


time 


ref. 






past experiments 








Heidelberg-Ms. 


76 Ge 


11 


ionization 


LNGS 


-2003 


[161 

i a 


Cuoricino 


iso Xe 


11 


bolometer 


LNGS 


-2008 


[36J 


NEMO-3 


100 Mo, 82 Se 


7,1 


track. +calorim. 


Modane 


-2011 


[37] 






current experiments 








EXO 


136 Xe 


175 


liquid TPC 


WIPP 


2011- 


1151 


Kamland-Zen 


136 Xe 


330 


liquid scintil. 


Kamioka 


2011- 


1 38 1 


GERDA-I/II 


76 Ge 


17/35 


ionization 


LNGS 


2011-/13 


[39J 


CANDLES 


48 Ca 


0.35 


scint. crystal 


Oto Cosmo 


2011- 


40 

L J 






funded experiments 








NEXT 


136 Xe 


100 


gas TPC 


Canfrac 


2014 


[4T1 


CuoreO / Cuore 


i30 Te 


10/200 


bolometer 


LNGS 


2012/14 


142 

i i 


Majorana Demo. 


76 Ge 


30 


ionization 


SUSEL 


2014 


03] 


SNO+ 


150 Nd 


44 


liquid scint. 


Sudbury 


2014 


m 






proposal, proto-typing 








SuperNEMO 


82 Se 


7/100-200 


track. +calorim. 


Modane 


2014/- 


m 


Cobra 


116 Cd 




solid TPC 


LNGS 




[46j 


Lucifer 


82 Se 




bolom.+scint. 


LNGS 




m\ 


DCBA/MTD 


i50 Nd 


32 


tracking 






[48j 


MOON 


82 Se, 100 Mo 


30-480 


track. +scint. 






[49j 


XMASS 


136 Xe 




liquid scint. 


Kamioka 




[50j 


AMoRE 


100 Mo 


100 


bolom.+scint. 


Yang Yang 




[51j 


Cd exp. 


116 Cd 




scint. 






[52j 



95 among the isotopes for a given method. The right hand plot of Fig [T] shows these ratios for the 

96 different models normalized to the decay rate of 130 Te. One sees that' 6 Ge is less favorable. The 

97 expected decays per kg vary between 20% and 50% of the rate of 130 Te. In other words: if all 

98 experimental parameters were the same then one would need a factor of ~3 more target mass in 

99 a 76 Ge experiment to have the same sensitivity. In reality this is not the case, i.e. the superior 

100 energy resolution of Ge detectors compensates this effect. 

101 5. Comparison of experiments 

102 The experiments searching for Ov/3/3 decay use a large variety of detection mechanisms and 

103 background reduction methods, see Tab. [2] The current status of almost all of them is described 

104 in these proceedings. Therefore a more detailed discussion is omitted here. Instead the key 

105 performance numbers are taken for a comparison of the sensitivities of some experiments. 

106 Since experiments use different isotopes a relative scaling factor for the different matrix 

107 elements and phase spaces has to be applied. This factor can be estimated using Fig. [TJ The 

108 values used here are /^(Ge) = 0.35, /^(Se) = 1.1, /^(Mo) = 1.6, /^(Te) = 1, /^(Xe) = 0.55 

109 and /A(Nd) = 2.2. 

If the number of background events is large, equation [8] can be used to estimate the 
experimental sensitivity. A relative figure-of-merit can then be defined as 

FoM =^-"-vS 

no One can call this the "ultimate" relative sensitivity of an experiment. Tab. [3] lists the 



(10) 



Table 3. Comparison of flgure-of-merits (FOM) for the case of large number of background 
events ("ultimate sensitivity"), /a is the scale factor for a given isotope taken from Fig.[T](right), 
and AE is the energy window which is taken here to be 1(2) full width half maximum for 
experiments with > 0.5% (< 0.5%) resolution. Note that the efficiency is reduced by 0.7 if 

AE = 1-FHWM. FOM is defined in the text. 

experiment mass Ja background AE efficiency enrichment FOM 





[kg] 




r cnt I 
1-kcV-kg-y-l 


[keV] 








Hd-Moscow 


11 


0.35 


0.12 


8 


0.8 


0.86 


0.8 


Cuoricino 


41 


1 


0.16 


12 


0.9 


0.27 


1.1 


NEMO-3 


6.9 


1.6 


0.002 


240 


0.18 


0.9 


1.0 


EXO 


175 


0.55 


0.004 


260 


0.33 


0.81 


1.9 


Kamland-Zen 


330 


0.55 


0.0002 


250 


0.5 


0.9 


20 


GERDA-I 


15 


0.35 


0.03 


10 


0.8 


0.86 


1.7 


GERDA-II 


30 


0.35 


0.001 


6 


0.8 


0.88 


17 


Major. -Dem. 


20 


0.35 


0.001 


6 


0.9 


0.9 


16 


CUORE 


750 


1 


0.01 


10 


0.9 


0.27 


21 


SNO+ 


800 


2.2 


0.0002 


230 


0.33 


0.056 


5.4 


NEXT 


100 


0.55 


0.0002 


25 


0.25 


0.9 


18 


SuperNEMO 


100 


1.1 


0.0002 


120 


0.3 


0.9 


19 


Lucifer 


100 


1.1 


0.001 


10 


0.9 


0.5 


50 



performance numbers and the figure-of- merit. For running (and past) experiments like EXO 
and GERDA-I the current achieved values are used which might improve with time while for 
the others the anticipated performance numbers are takenj^] 

Alternatively, the (relative) sensitivity vs. time can be estimated from equation [7] by 



f0u f A -e- V -M-t 
Jl / 2> V(B-AE-M-t) 



lu Here ^(Abkg) is the "average" 90%C.L. upper limit of the number of signal events for Abkg 

us background events calculated according to the method discussed in [29] . The result is shown in 

116 Fig [2} Here all experiments are assumed to start at time 0. 

117 A few comments should be made concerning the interpretation of Tab. [3] and Fig. [2j 

us • If one takes the spread of the data points in Fig. [T] the factor /a has a ~ 20% uncertainty. 

119 • The 2^/3/3 background is irreducible and can only be avoided with an energy resolution 

120 a < 1 — 2% at Qpp. This requirement depends of course strongly on Ty 2 which varies by 

121 a factor of 300 for the isotopes considered. For some experiments this background is not 

122 fully taken into account for the background index. 

123 • All sensitivities given are the scales for 0^(3(3 discovery. To get relative sensitivities for mpp 

124 the square root has to be taken. 

125 • Of the running experiments, Kamland-Zen should have the largest potential. This is 

126 impressive if one takes into account that it was not specially built for this physics. 

127 • Germanium experiments can be very competitive despite the fact that the phase space 

128 factor is so small. Especially if a positive signal will be claimed, a narrow peak at Qpp will 

129 be more convincing than a broad shoulder. 

2 A fiducial volume cut will reduce the active mass. Depending on whether the background index is normalized 
to the total mass or to the fiducial mass, the efficiency e has to go under the square root or not. The meaning of 
B is not always clearly defined in the literature. Here the normalization to the total mass is assumed. 
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Figure 2. Relative experimental sensitivity for 0^/3/3 life time limit versus running time. 



130 • The Lucifer approach with 100 kg is very competitive even in comparison to a ton scale Xe 

131 experiment like Kamland-Zen or NEXT. 

132 • Systematic effects like the precision of the energy resolution or the background shape are 

133 not taken into account. 

134 In case the neutrino masses are ordered in the inverted mass hierarchy, a lower bound of 

135 about 15 meV for mpn can be calculated using the current parameters from neutrino oscillation 

136 experiments. For Ge this corresponds to half lives of 5 — 20 • 10 27 years. These values should 

137 be compared to the expected sensitivity of GERDA-II or the Majorana Demonstrator of about 

138 1.5-10 26 y. This demonstrates that exploring the entire mass band of the inverted hierarchy is a 

139 long term enterprise. With the numbers in Tab. [3] and a mass of 1000 kg, the required time for 
wo 5 • 10 27 y is 13 years while a Lucifer like experiment would need to run for half the time. 

141 6. Summary 

142 Neutrinoless double beta decay is the best experimentally accessible method to test whether 

143 neutrinos are Majorana particles. This decay violates lepton number and is therefore on equal 

144 footing to proton decay searches. The motivation for several large efforts in this field is therefore 

145 obvious. 

146 For a long time, the Heidelberg-Moscow experiment has dominated the field and its claim of 

147 a Oz^/3/3 signal has not been scrutinized since 2001. In 2011, EXO, Kamland-Zen, CANDLES and 

148 GERDA-I started to take data. All but CANDLES are more sensitive than Heidelberg-Moscow 

149 and especially Kamland-Zen is expected to answer this question in the next 12 months. EXO has 

150 already reported a first time measurement of T 2 ^ 2 ( 136 Xe) = 2.11 ± 0.04(stat) ± 0.21(syst) • 10 21 y 

151 which is considerably lower than previous limits |15| . 

152 Beyond this next step, experiments want to explore the map region for the inverted neutrino 

153 mass hierarchy. This will eventually require ton scale experiments. Which of the proposed 

154 solutions will be built is open at the moment. 
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